Incommensurate magnetic fluctuations and Fermi surface topology in LiFeAs 
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Using the angle-resolved photoemission spectroscopy (ARPES) data accumulated over the whole 
Brillouin zone (BZ) in LiFeAs we analyze the itinerant component of the dynamic spin susceptibility 
in this system in the normal and superconducting state. We identify the origin of the incommensu- 
rate magnetic inelastic neutron scattering (INS) intensity as scattering between the electron pock- 
ets, centered around the (7r,7r) point of the BZ and the large two-dimensional hole pocket, centered 
around the F-point of the BZ. As the magnitude of the superconducting gap within the large hole 
pocket is relatively small and angle dependent, we interpret the INS data in the superconducting 
state as a renormalization of the particle-hole continuum rather than a true spin exciton. Our 
comparison indicates that the INS data can be reasonably well described by both the sign changing 
symmetry of the superconducting gap between electron and hole pockets as well as sign preserving 
gap, depending on the assumptions made for the fermionic damping. 

PACS numbers: 74.70.Xa, 74.20.Fg, 75.10.Lp, 75.30.Fv 



The relation between unconventional superconductiv- 
ity and magnetism is one of the most interesting topics 
in condensed-matter physics. For example, in most of 
the iron-based superconductors superconductivity occurs 
in close vicinity to an antiferromagnetic (AF) state [1- 
3]. Moreover, superconductivity emerges when antiferro- 
magnetic order in parent compounds is suppressed, either 
by electron/hole doping or disorder. In addition, short- 
range AF spin excitations are still present in the normal 
state of the doped systems and also become resonant in 
the superconducting state at energies below twice the 
superconducting gap magnitude, 2Ao[4]. This resonant 
enhancement is believed to be a signature of a certain 
phase structure of the superconducting gap as the para- 
magnetic spin response of the Bogolyubov quasiparticles 
at the antiferromagnetic wave vector Qaf sensitive to 
the anomalous coherence factor 1 — Ak+q^^ ^ Once 
the superconducting gap at parts of the Fermi surface, 
connected by Qaf, changes sign, the spin response ac- 
quires an additional enhancement at 1^ < 2Ao, which 
is a hallmark of unconventional superconductivity. The 
observation of the spin resonance in many iron-based su- 
perconductors provides strong evidence for the so-called 
-wave symmetry of the superconducting gap, where 
the gap structure changes sign between electron and hole 
pockets [5-7]. Note that this does not exclude the gap on 
each pocket to have a strong angular variation and even 
accidental nodal lines, allowed by Aig symmetry [3]. The 
angular variation of the gap, measured in ARPES [8], is 
inconsistent with idealized lattice version of 5^ , but can 
be modeled by taking into account interaction effects. 



While the behavior, described above, is observed in 
the majority of the iron-based superconductors, there 
are some notable exceptions. Perhaps the most inter- 
esting one is the stoichiometric LiFeAs, which supercon- 
ducts at Tc=17K without any doping[9-ll]. In addi- 
tion, LiFeAs shows neither static AF ordering nor nest- 
ing between electron and hole bands at Q^i?[12]. Several 
neutron scattering experiments were performed recently 
in LiFeAs [13-15], including only one study [15] on super- 
conducting single crystals, where magnetic intensity at 
an incommensurate momentum close to (tt, tt) was ob- 
served. Its renormalization across Tc was found to be 
too weak to draw a definite conclusion about the phase 
structure of the superconducting gap. Furthermore, some 
controversy on the phase structure of the order parame- 
ter arises in the analysis of the quasiparticle interference 
in the superconducting state of LiFeAs[16, 17]. 

In this paper we use the ARPES data for LiFeAs, 
Ill-type pnictide superconductor, which is known to be 
free from surface effects [19] to strengthen a connection to 
the INS response. This is particularly important, given 
the controversy on the Fermi surface topology in this 
system. [12, 18] We employ an effective tight-binding fit to 
the high-quality LiFeAs photoemission data in order to 
compute the spin response within random phase approxi- 
mation (RPA). We believe this procedure is only possible 
at present in 111 systems as availability of the requi- 
site data is most complete here. A comparison with INS 
data shows that the incommensurate magnetic scatter- 
ing intensity arises due to scattering between the electron 
pockets, centered around the (7r,7r) point of the BZ, and 
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the large two-dimensional hole pocket, centered around 
the F-point of the BZ. We also find that the renormal- 
ization of the neutron intensity upon opening of the su- 
perconducting gap is relatively weak, consistent with the 
INS experiments. 

ARPES and tight-binding fit. Owing to its ability to re- 
solve both momentum and energy of the electronic states, 
modern photoemission can be used to map out a com- 
plete low-energy electronic structure of a layered com- 
pound, like LiFeAs. From such a comprehensive data 
set one may extract the dispersion of quasiparticles at 
any momentum. This can be used to calculate numerous 
properties like heat capacity, plasma frequency or the 
Hall coefficient [20, 21]. However to make this possible 
quasiparticle dispersions have to be conveniently param- 
eterized. One way to do this is via a tight-binding fit. 
Indeed, tight-binding models including up to 10 bands, 
have been developed to fit the LDA band structure of the 
iron-pnictide superconductors [22, 23]. However from a 
practical perspective it is more favorable to use an ef- 
fective tight-binding model, separately describing disper- 
sions of each band that crosses the Fermi level [24]. In 
case of the square lattice with a tetragonal symmetry 
the quasiparticle dispersion can be fit by the following 
formula: 
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FIG. 1. (color online) (a) Experimental Fermi surface map 
with fitting contours superimposed over it. (b) Three gen- 
eral energy-momentum cuts passing through the electron 
pocket centered at (— vr/a; — 7r/5) and hole pockets at (0; 0) 
and (0;-27r/6). 



where am,n is di N x N matrix of effective tight-binding 
coefficients, and (l)m,n{kx^ky) are base functions 
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The a matrices are chosen to provide the best fit for the 
form of the Fermi surface pockets and band velocities at 
the Fermi level. For the two hole pockets centered at the 
F point and for the two electron pockets located at the 
corners of the BZ, the parameters are as follows (in eV): 
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To demonstrate to which extent this simple model is 
able to capture the dispersion of low energy bands, in 
Fig. 1 we plot typical experimental FS maps[12, 19, 20] 
with the fitted band dispersion. The lower panel also 
contains several general energy-momentum cuts, allow- 
ing one to compare the experimental and model disper- 
sions for energies close to the Fermi level. 

Spin response. In the following we proceed with the 
calculations of the magnetic INS spin response. In the 
magnetically-disordered state transverse and longitudi- 
nal components of the spin susceptibility are identical, 
and we focus below on the transverse part. The spin 
response is computed within RPA. Then, the transverse 
components of the full spin susceptibility x*'"^ are related 
to transverse components of the bare susceptibility Xo"^ 
as 
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Although the electron bands are known to have notice- 
able kz dispersion we use averaged parameters for the 
sake of simplicity, thus remaining within a 2D structure. 



where i and j are band indices. 

Summation over repeated band indices is implied and 
^ are matrix elements of the interactions. The solution 
of Eq.(3) in matrix form is straightforward: x = Xo(l — 
uxq)~^ . The components of the bare spin susceptibility 
Xo = Xo"^(Q7i^m) are given by usual combinations of 
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normal and anomalous Green's functions 
T 

"27V 
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where Gl.{iuJn) = 

At 



, M^IYHK)' and F^(iw„) = 
For the superconducting gap func- 



tion we assume the form, obtained in the ARPES 
experiments [8]. In particular, the superconducting gaps 
on the inner and on the outer hole pockets were found 
to amount to /\^^^^^^ = 6meV, and A^°^*^^ = 3.4 + 
0.5 (cos 4^ + 0.13 cos 80 — 0.2 cos 12^) (in meV), respec- 
tively. Here (j) is the angle counted on the hole Fermi sur- 
face. For the two electron pockets the gaps were found to 
be similar in the form A^*'"^^'^ = A^^"*^'^ = 3.4 + 0.5 cos 6 
(in meV) where is the angle on the electron Fermi sur- 
face pockets. Note that a similar angular variation of the 
superconducting gap on the outer hole pocket was ex- 
tracted in Scanning Tunneling Microscopy (STM) though 
with a smaller gap magnitude. [17]. As ARPES is not sen- 
sitive to the phase difference of the gap between electron 
and hole pockets we consider two possibilities, namely, 
5^ -symmetry of the superconducting gap, where the 
phase of the superconducting gap changes sign between 
electron and hole pockets and s++ where, despite higher 
harmonics, the gaps on the electron and hole pockets re- 
mains always positive. 

In our numerical calculations we keep all terms in the 
matrix equation for the full susceptibility. The interact- 
ing part of the Hamiltonian contains four-fermion inter- 
actions with small momentum transfer as well as mo- 
mentum transfers around (7r,7r). They include the inter- 
actions between electron and hole bands with momentum 
transfer around (tt, tt) as well as interactions with small 
momentum transfer within or between hole pockets and, 
similarly, for the electron pockets. For simplicity, we ap- 
proximate all interactions as angle-independent, i.e., we 
neglect the angle dependence introduced by dressing the 
interactions by coherence factors associated with the hy- 
bridization of Fe (i— orbitals. These coherent factors do 
play a role in the angular variation of the superconduct- 
ing gap [25], but do not substantially modify the positions 
of the spin resonance [26]. For better convergence of the 
numerical series we add a small damping of F = 3 meV 
to the fermionic dispersion in the normal state. This 
value is consistent with values, extracted from ARPES 
experiments [12]. We also analyze the influence of the 
fermionic damping on the spin excitations in the super- 
conducting state below in more detail. 

We start by looking on the bare susceptibility as this 
quantity directly follows from the fermiology, measured 
by ARPES. In particular, in Fig. 2 and Fig. 3 (a) we show 
its imaginary part as q and Vt maps in the first BZ. As ex- 
pected, the scattering momenta associated with 2kF in- 



traband processes resemble the original Fermi surfaces in 
circular-like structures. Furthermore, Fig. 2 allows for a 
straightforward identification of the character of the scat- 
tering. The intraband scattering and interband scatter- 
ing between two electron or two hole bands are centered 
around (0, 0), while the interband scattering between the 
bands of different character are centered around (±7r, ^tt) 
momentum. By comparing the diameters of the intra- 
band driven scattering circles with the approximate 2kF 
values of the corresponding Fermi surfaces one identifies 



immediately 2kp'^*'''^ as well as 2k^ 
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shown 



by the arrows. In addition, the bright spot around the 
q = refers to the scattering within the small inner hole 
pocket as well as the scattering between the inner and 
the outer electron pockets. Further interband scattering 
processes include the scattering between two hole pock- 
ets, denoted by qhh and most importantly, the scatter- 
ing between the electron pockets and outer hole pocket, 
shown by Q^. Note that the scattering between the in- 
ner hole pocket and the electron pockets occurs also at 
a similar momentum, but its intensity is much smaller 
already in the bare susceptibility due to limited phase 
space available for scattering. Within the RPA the in- 
tensity of these excitations is further suppressed as com- 
pared to the scattering between the outer hole pocket 
and electron pockets. Therefore, we can safely conclude 
that the excitations at arise due to scattering between 
the large outer hole pocket and the two electron pockets. 
Note, our analysis does not include the matrix elements, 
originating from the transformations from the orbital ba- 
sis to the band ones. Its inclusion usually strengthens the 
transverse scattering, q^^ = (g, 27r — g), over the longitu- 
dinal ones q^^ = {q^q)[27]. This would make the inten- 
sity, shown in Fig. 2, look more anisotropic, but it will not 
change the position of the peaks. Note that our calcu- 
lations within a three-orbital model are consistent with 
this observation. 

These scattering wave vectors are also visible along the 
transverse direction, qtr = (^, 27r — g), shown in Fig. 3(a). 
The scattering within the small inner hole pocket at small 
momentum q as well as the scattering between the elec- 
tron pockets and outer large hole pocket at the wave vec- 
tor are most pronounced. One can further identify the 
intraband scattering within the electron bands at 2k p. 
Note also that on the energy scale from to 15 meV the 
dispersion of these excitations is almost vertical, which 
is caused by the relatively large Fermi velocity of the in- 
volved bands. By comparison with the INS we find that 
the scattering momentum ^ (0.86, 1.14)7r, associated 
with the scattering between the outer hole pocket and the 
two electron pockets, matches precisely with the experi- 
mentally observed incommensurate momentum. As our 
band structure results from the fit to the ARPES band 
structure, we conclude that the incommensurate momen- 
tum seen in the INS refers to the scattering between the 
electron and the outer hole Fermi surfaces. This can be 
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FIG. 2. (color online) Calculated imaginary part of the bare 
spin susceptibility in the normal state of LiFeAs as a function 
of the momentum in the first BZ at /iQ=5meV. 
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FIG. 3. (color online) Calculated imaginary part of the bare 
(a) and the RPA (b) spin response as a function of the trans- 
verse momentum and frequency in the normal state of LiFeAs. 

further supported by the fact that the incommensurate 
magnetic excitations, found in INS, do not indicate a 
strong z-dispersion. We recall that the small inner hole 
pocket around the F— point of the BZ has strongly three- 
dimensional character which we ignored at present. How- 
ever, if taken into account, it should produce a strong 
dispersion of the incommensurate magnetic excitations 
along the Qz momentum, which is not the case. There- 
fore, the scattering between the outer hole pocket and 
the two electron pockets is, most likely, responsible for 
the INS intensity at the wavevector ^ (0.86, 1.14)7r, 
which differs from the proposal made in Ref.[28] where 
these incommensurate peaks were attributed to the scat- 
tering between the small inner hole pocket and two elec- 
tron pockets. 

To proceed further, we compute the total RPA suscep- 
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FIG. 4. (color online) Calculated momentum dependence of 
the the imaginary part of the RPA spin response at HQ. — 
5meV in the normal and suerconducting -wave states, re- 
spectively. The symbols refer to the experimental data, taken 
from Ref.[15]. 

tibility by including the interactions. Note, however, that 
part of these interaction parameters contributes already 
to the renormalization of the bands, which was used to 
obtain the tight-binding model. As the interaction val- 
ues are not fully known we took them to be the same, 

i.e. U = Uinterband = Uinterband = O.TSOf^^^l^ (1 , 2) CXCCpt 

for the scattering between the electron pockets, which 
we consider to be small Uee = O.lu. The magnitudes 
of the interactions were also chosen such that the sys- 
tem remains in the paramagnetic phase. Fig. 3(b) shows 
the results for the Iirxrpa^ displayed as qtr and E = HQ 
map. Observe that in comparison to the bare susceptibil- 
ity the incommensurate excitations due to the scattering 
between electron and hole bands are enhanced. At the 
same time, we find that the excitations at small q are 
much less intense as compared to the bare xo- To under- 
stand the origin of its suppression at finite frequencies 
recall that the real part of the bare intraband suscepti- 
bility falls off as 1/Q which indicates that within RPA 
there is no source for the enhancement of these small q 
excitations. This explains why the total susceptibility 
shows stronger enhancement only for the wavevector 
and not for the q ~ momentum. 

A comparison to the experimental INS data in the nor- 
mal state is shown in Fig. 4 where we display the imagi- 
nary part of the total RPA susceptibility for HQ = 5meV 
as a function of the transverse momentum. The pro- 
nounced peaks at due to the scattering between the 
outer hole and two electron pockets agree with those 
found experiment ally [15]. Nevertheless, one should men- 
tion that the structure of the peaks is symmetric in the 
calculations, while in INS there is an additional shoul- 
der for larger qtr- A weak 2:-dispersion or some other 
scattering paths may cause this behavior. 

In the next step, we move to the superconducting state 
and compute the spin excitations for various symmetries 
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FIG. 5. Calculated frequency dependence of the imaginary 
part of the RPA spin response at the wavevector in the nor- 
mal state and the superconducting state for two different sym- 
metries of the superconducting gap. Difference curves for the 
s~^~^ and -wave scenarios refer to the constant fermionic 
damping of 3 meV and the one frequency dependent, as de- 
scribed in the text. The inset represent the experimental data, 
taken from Ref.[15]. Observe that the upturn behavior of the 
experimental curves at energies smaller than 1.5 meV is due 
to elastic scattering contribution. 

of the superconducting order parameters. The most in- 
teresting question is whether any information can be ex- 
tracted about the phase structure of the gap with respect 
to the relative phase difference between electron and hole 
pockets as well as between inner and outer hole pockets. 
We remind that in contrast to the angular dependence 
of the gap, the relative phase structure cannot be di- 
rectly probed by ARPES. In particular, we considered 
two different situations. The first one, which we name 
, refers to the phase of the superconducting gap on 
the hole pockets being opposite to the phase of the gap 
on the electron pockets. In the other case, the so-called 
5++, the overall phase of the order parameter is the same 
for the electron and the hole pockets. 

In Fig. 4 we show the behavior of spin excitations for 
the -wave symmetry together with the normal state 
results and the experimental data for h^} = 5meV. Ob- 
serve that the renormalization of the spin excitations in 
the 5^ —wave channel is present but relatively weak in 
the sense that the excitations are only slightly enhanced 
with respect to the normal state. The reason for this 
moderate renormalization is the relatively strong angu- 
lar variation of the superconducting gap on the outer hole 
pocket and on the electron pockets [29]. The angular- 
dependent gap washes out the strong enhancement of 
Imxo at 2Ao which is a prerequisite for the sharp reso- 
nance. In LiFeAs the enhancement of Imxo in the super- 
conducting state with respect to its normal state value 
is more gradual than in the other FeAs superconduc- 
tors. In addition the gaps on the electron and on the 
outer hole pockets contributing mostly at this Qi are rel- 



atively small. Taken together with the constant value 
of the fermionic damping of the order of ~ 3meV, these 
factors render the enhancement of the spin excitations in 
the -wave scenario relatively weak. In other words 
in LiFeAs there is no true spin resonance below the con- 
tinuum of particle-hole excitations but an enhancement 
of the continuum itself due to coherence factors associ- 
ated with the phase structure of the -wave supercon- 
ducting gap. This also makes it difficult to distinguish 
the other scenarios. In particular, although an -wave 
gap structure produces spin excitations which agree with 
available experimental data, the results with an 5++- 
wave symmetry cannot be ruled out based on the mea- 
surements at a given frequency. The difference between 
two symmetries of the order parameters becomes more 
apparent in Fig. 5 where we plot Im Xrpa = Xli j X*'"^ 
the superconducting and in the normal state, as a func- 
tion of frequency (HQ) at the wavevector Q^. Notice first 
that the normal state curve shows a characteristic sin- 
gle relaxor form of the overdamped paramagnons, cen- 
tered around 7.5 meV. It agrees qualitatively with the 
INS data, see Fig.3(b) in Ref.[15]. 

In the superconducting state the results for the s++ 
and the s^ -wave superconducting gaps depend sensi- 
tively upon the assumption made for the fermionic damp- 
ing. For constant damping the is the only sym- 
metry which qualitatively agrees with the experimental 
data of Ref.[15], shown in the inset of Fig. 5. In par- 
ticular, one finds that the intensity of the spin excita- 
tions is suppressed with respect to its normal state val- 
ues up to energies of about 4.2 meV and is then slightly 
enhanced for higher energies. However, as mentioned 
above this enhancement is not a true exciton but rather 
an enhancement of the particle-hole continuum due to 
the sign change of the superconducting gap. At the 
same time, for the s++-wave symmetry the difference 
between the superconducting and normal states remains 
always negative for energies up to 30 meV for constant 
damping. Here, the spin excitations are suppressed as 
the superconducting gap does not change sign at this 
particular Q^. We note, however, that manipulation of 
the fermionic damping improves the situation. Follow- 
ing Ref.[30] we modeled the fermionic damping in the 
superconducting state as F ~ for < HQ < SAave^ 
F = Tns for an HQ > AAave^ and increasing linearly for 
SAave < < 4:Aave- We Varied the value for Aave be- 
tween 3 and 6 meV to find the best-case scenario for the 
5++-wave symmetry regarding the enhancement in the 
superconducting state. The result for Aave ^ 5meV is 
shown in Fig. 5. For the chosen damping parameters, the 
s++-wave symmetry exhibits scattering enhancement in 
the superconducting state while a true spin exciton at the 
energy of about 6meV appears for the s^ -wave symme- 
try. However, the enhancement for s++ occurs at higher 
energies than in the experiment which clearly indicates 
intensity gain in the superconducting phase for energies 
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between 5 and 10 meV. The agreement with s++-wave 
symmetry can be improved by reducing the Aave value, 
but the magnitude of the enhancement with respect to 
the normal state then also becomes smaller. At this 
point, we can only conclude that experimental results do 
not show a true spin resonance mode, frequently taken 
as characteristic for the s^ -superconductor, but a redis- 
tribution of the particle-hole continuum in the presence 
of superconductivity. This renders the definite conclu- 
sion on the phase structure of the superconducting gap 
rather difficult. 

To conclude, we employ the ARPES data for LiFeAs 
to make a connection to the INS response. A compar- 
ison to the INS data shows that the incommensurate 
magnetic scattering intensity originates from the scatter- 
ing between the electron pockets, centered around the 
(tt, tt) point of the BZ and the large two-dimensional 
hole pocket, centered around the F-point of the BZ. This 
points towards an internal consistency between the FS 
topology, measured by ARPES, and INS results. We also 
find that the renormalization of the neutron intensity in 
LiFeAs in the superconducting state can be understood 
in terms of a rearrangement of the particle-hole contin- 
uum which is rather weak for any phase structure of the 
superconducting gap between electron and hole pockets. 
Further studies are necessary to elucidate the nature of 
the Cooper-pairing in this compound. 

We thanks J. van den Brink, A.V. Chubukov, and 
P.J. Hirschfeld for fruitful discussions. IE acknowl- 
edges financial support of the Merkur Foundation, Ger- 
man Academic Exchange Service (DAAD PPP USA No. 
50750339). JK acknowledges support from the Studien- 
stiftung des deutschen Volkes, IMPRS Dynamical Pro- 
cesses in Atoms, Molecules and Solids, and DFG within 
GRK1621. VZ and SB acknowledge support by DFG 
grants ZA 654/1-1 and BO 1912/3-1 (SPP1458). NQ and 
MB were supported by the DFG through SFB 608. The 
work is supported by DFG grants BE1749/13 (SPP1458) 
and the Graduiertenkolleg GRK1621 and by the Dresden 
Platform for Superconductivity and Magnetism. 



* jknolle@pks.mpg.de 
^ ieremin@tp3.rub.de 
[1] Y. Kamihara, T. Watanabe, M. Hirano and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] P. J. Hirschfeld, M. M. Korshunov, 1. 1. Mazin, Rep. Prog. 

Phys. 74, 124508 (2011). 
[3] A. V. Chubukov, Annu. Rev. Condens. Matter Phys. 3, 
13 (2012). 

[4] M. D. Lumsden and A. D. Christianson, J. Phys.: Con- 
dens. Matter 22, 203203 (2010). 

[5] K. Kuroki et al., Phys. Rev. Lett. 101, 087004 (2008). 

[6] I.I. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 

[7] A. V. Chubukov, D. V. Efremov, and 1. Eremin, Phys. 



Rev. B 78 134512 (2008). 
[8] S.V. Borisenko, V. B. Zabolotnyy, A.A. Kordyuk, D.V. 

Evtushinsky, T.K. Kim, I.V. Morozov, R. Follath, and 

B. Biichner, Symmetry 4(1), 251 (2012) 
[9] X.C. Wang, Q.Q. Liu, Y.X. Lv, W.B. Gao, L.X. Yang, 

R.C. Yu, F.Y. Li, and C.Q. Jin, Solid State Commun. 

148, 538 (2008). 
[10] C. W. Chu, F. Chen, M. Gooch, A.M. Guloy, B. Lorenz, 

B. Lv, K. Sasmal, Z.J. Tang, J.H. Tapp, and Y.Y. Xue, 

Physica C 469, 326 (2009). 
[11] I. Morozov, A. Boltalin, O. Volkova, A. Vasiliev, O. 

Kataeva, U. Stockert, M. Abdel-Hafiez, D. Bombor, A. 

Bachmann, L. Harnagea, M. Fuchs, H.-J. Grafe, G. Behr, 

R. Klingeler, S. Borisenko, C. Hess, S. Wurmehl, and B. 

Biichner. Cryst. Growth and Design 10, 4428 (2010). 
[12] S.V. Borisenko, V.B. Zabolotnyy, D.V. Evtushinsky, 

T.K. Kim, I.V. Morozov, A.N. Yaresko, A.A. Kordyuk, 

G. Behr, A. Vasiliev, R. Follath, and B. Biichner, Phy. 
Rev. Lett. 105, 067002 (2010). 

[13] A.E. Taylor, M.J. Pitcher, R.A. Ewings, T.G. Perring, 
S.J. Clarke, and A.T. Boothroyd, Phys. Rev. B 83, 
220514 (R) (2011). 

[14] M. Wang, X.C. Wang, D.L. Abernathy, L.W. Harriger, 

H. Q. Luo, Y. Zhao, J.W. Lynn, Q.Q. Liu, C.Q. Jin, C. 
Fang, J. Hu, and P. Dai, Phy. Rev. B 83, 220515(R) 
(2011). 

[15] N. Qureshi, P. Steffens, Y. Drees, A.C. Komarek, D. 
Lamago, Y. Sidis, L. Harnagea, H.-J. Grafe, S. Wurmehl, 
B. Biichner, and M. Braden, Phys. Rev. Lett. 108, 
117001 (2012). 

[16] T. Hanke, S. Sykora, R. Schlegel, D. Baumann, L. Har- 
nagea, S. Wurmehl, M. Daghofer, B. Biichner, J. van den 
Brink, and C. Hess, Phys. Rev. Lett. 108, 127001 (2012). 

[17] M.P Allan, A.W. Rost, A.P Mackenzie, Y. Xie, J.C. 
Davis, K. Kihou, C.H. Lee, A. lyo, H. Eisaki, and T.-M. 
Chuang, Science 336, 563 (2012). 

[18] C. Putzke, A.I. Coldea, I. Guillamon, D. Vignolles, A. 
McCollam, D. LeBoeuf, M.D. Watson, LI. Mazin, S. 
Kasahara, T. Terashima, T. Shibauchi, Y. Matsuda, and 

A. Carrington, Phys. Rev. Lett. 108, 047002 (2012). 
[19] A. Lankau, K. Koepernik, S. Borisenko, V. Zabolotnyy, 

B. Biichner, J. van den Brink, and H. Eschrig, Phys. Rev. 
B 82, 184518 (2010) 

[20] U. Stockert, M. Abdel-Hafiez, D.V. Evtushinsky, V.B. 

Zabolotnyy, A.U.B. Wolter, S. Wurmehl, I. Morozov, R. 

Klingeler, S.V. Borisenko, and B. Biichner, Phys. Rev. B 

83, 224512 (2011) 
[21] D.V. Evtushinsky, A.A. Kordyuk, V.B. Zabolotnyy, D.S. 

Inosov, B. Biichner, H. Berger, L. Patthey, R. Follath, 

and S.V. Borisenko, Phys. Rev. Lett. 100, 236402 (2008) 

[22] H. Eschrig and K. Koepernik Phys. Rev. B 80, 104503 
(2009) 

[23] S. Graser, A.F. Kemper, T.A. Maier, H.-P. Cheng, P.J. 
Hirschfeld, and D.J. Scalapino, Phys. Rev. B 81, 214503 
(2010) 

[24] D. S. Inosov, V.B. Zabolotnyy, D.V. Evtushinsky, A.A. 

Kordyuk, and B. Biichner, R. Follath, H. Berger, and 

S.V. Borisenko, New J. Phys. 10, 125027 (2008). 
[25] S. Maiti, M.M. Korshunov, T.A. Maier, P.J. Hirschfeld, 

A.V. Chubukov, Phys. Rev. B 84, 224505 (2011) 
[26] S. Graser, T. A. Maier, P. J. Hirschfeld, and D. J. 

Scalapino, New J. Phys. 11, 025016 (2009). 
[27] J.T. Park, D.S. Inosov, A. Yaresko, S. Graser, D.L. Sun, 



7 



Ph. Bourges, Y. Sidis, Y. Li, J.-H. Kim, D. Haug, A. 
Ivanov, K. Hradil, A. Schneidewind, P. Link, E. Faul- 
haber, L Glavatskyy, C. T. Lin, B. Keimer, and V. 
Hinkov, Phys. Rev. B 82, 134503 (2010). 
[28] M. Wang, M. Wang, H. Miao, S.V. Carr, D.L. Abernathy, 
M.B. Stone, X.C. Wang, L. Xing, C.Q. Jin, X. Zhang, 



J. Hu, T. Xiang, H. Ding, and P. Dai, arXiv: 1208.0909 
(unpubhshed) 

[29] S. Maiti, J. Knolle, L Eremin, and A.V. Chubukov, Phys. 

Rev. B 84, 144524 (2011). 
[30] S. Onari, and H. Kontani, Phys. Rev. B 84, 144518 

(2011) . 



